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Subwavelength Transportation of Light with Atomic Resonances 
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We propose and investigate a new type of optical waveguide made by an array of atoms without 
involving conventional Bragg scattering or total internal reflection. A finite chain of atoms collec¬ 
tively coupled through their intrinsic resonance supports a propagating mode with minimal radiative 
loss when the array spacing a is around 0.6Ao/27r where Ao is the wavelength of the nearly resonant 
optical transition. We find that the transportation is robust with respect to position fluctuation 
and remains possible when the atoms are placed on a circle. Our result paves the way to implement 
the subwavelength transportation of light in integrated optical circuits with cold atoms. 


There is much recent development of integrated opti¬ 
cal circuits localising and interfacing cold atoms in which 
quantum networks, quantum many-body physics, and 
quantum optics can be explored HHa. This new in¬ 
terdisciplinary toolbox ranging from atomic physics and 
quantum optics to nanophotonics opens the new door to 
investigate a wide range of remarkable phenomena such 
as the coherent transportation of photons emitted from 
an atom 016 ] and atomic mirrors [7]. In this Letter, 
we examine a waveguide made with an array of atoms 
along linear and circular geometries. Nanoscale waveg¬ 
uides that confine, guide and manipulate the electromag¬ 
netic energy on subwavelength scales have been actively 
studied as optoelectronic circuits ISHIQ]. The plasmon- 
based waveguides such as those made with nanowires or 
metallic particles dlHIS] have very excellent lateral con¬ 
finement but the propagation distance is considerably 
limited. This originates from the intrinsic loss of the 
metal which hampers the implementation of nanoscale 
metal plasmon waveguides. This loss is reduced when di¬ 
electric scatterers are used m via so-called “geometric 
resonances”. For scatterers with size R, however, the ge¬ 
ometric resonance strictly requires the lattice spacing a 
satisfying 2R < a < X/2 where A is the wavelength of 
the guided light. Thus, the finite size R of the scatterer 
considerably limits the flexibility of the waveguide made 
by scatterers. It has not been possible to demonstrate 
waveguides with arbitrary curvature from rigorous first 
principles numerical calculations. Nor is the waveguid- 
ing properties of simple arrays of scatterers robust under 
positional fluctuations. 

Here, we consider a waveguide made with an array 
of atoms where there is no size restriction. Instead of 
the “geometric resonance” with dielectric scatterers, we 
employ intrinsic “atomic resonance” to support a propa¬ 
gating mode in flexible sub-wavelength geometries in this 
Letter. With appropriate lattice spacing a ^ 0 . 6 Ao/ 27 r 
of the atomic chain where Aq denotes the wavelength of 
nearly resonant light, the external electric field exciting 
the first few atoms can propagate along the finite size of 
the chain with minimal radiation loss. We find that the 
transportation is robust with respect to position fluctua¬ 
tion of atoms and is possible when the atoms are placed 


on a circle. Recent development of plasmonic nanosphere 
[18] or two-dimensional photonic crystals [19] localising 
atoms should allow one to prepare an array of atoms with 
tunable lattice spacing. 

Through the excitation in an atom at one end of the 
chain, energy is transfered to all adjacent atoms in turn 
through the resonant excitations. For the chain of atoms, 
resonant excitations in different atoms are coupled by 
electromagnetic interaction and form collective optical 
modes. We consider the transportation of this excitation 
along the atomic chain is through an intrinsic p-wave 
resonance. Then, we investigate how excitation of the 
first atom by an external electric field E can travel down 
the chain. 

We expand the field near each atom locally in terms 
of vector spherical harmonics and examine the scattering 
phase shift of electric fields from each atom as is described 
in Ref. [20]. To identify the waveguide mode along the 
chain of atoms, we employ the multiple scattering theory 
describing electromagnetic interactions between the dif¬ 
ferent atoms [21]. Here, the scattered electric field from a 
single atom is related to the sum of the external electric 
field and the scattered fields from all the other atoms. 
As we see below, the waveguide mode is indicated by the 
undiminished propagation of the electromagnetic wave 
down the chain after the first few atoms are excited. 

In the multiple scattering approach EB, the total in¬ 
coming electric field driving the i^^ atom is the 

sum of the external electric field E?^^ and the scat¬ 
tered fields from all the other atoms: '^j^i ie 

= E?®* + illustrated in Fig. 0 This 

incoming field is scattered by the atom and pro^ces a 
scattered field Ethat is equal to the t matrix times the 
incoming field: E-^^ = tE-^\ The field E-^^^is equal to 
the photon propagator times the scattered field from 

site j: E-^^^ = GijEj^\ Here, the t matrix is diagonal 
due to the rotational invariance of the Hamiltonian [20] , 
and given by t = tan^/(tan(5 + i) where the scattering 
phase shift satisfies tan (5 = 1 /{uj — ujq P iV) for the natu¬ 
ral linewidth F of the optical transition. To summarise, 
we get the set of linear equations 
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FIG. 1: (color online) Our modelled system. An array of 
two-level atoms responding to the external electromagnetic 
field at the frequency of u. The atom is drived by 
both external fields and scattered fields from atoms 
where j ^ i. 


Ef ^ GijEf ^ + Ef *] (1) 


In momentum space, the multiple scattering equation 


E(3)()!c) 


1 — t(k)G(k) 


( 2 ) 


gives the resonance condition t(k)G(k) = 1. 

Linear Chain To illustrate the subwavelength trans¬ 
port of photons, we first consider a one-dimensional array 
of 30-320 atoms with a p-wave resonance at an angular 
frequency uoq and separated by lattice spacings a. The 
distribution of the electric field as a function of the posi¬ 
tion of the atoms is numerically calculated when the first 
few atoms are excited by the external electric field. We 
find that there can be a large response at the end of the 
chain but only when the external field is perpendicular 
to the line of atoms. We attribute this to the nature of 
p-wave resonance of atoms. In addition, we investigate 
the dependence on the way the system is excited. The 
propagated waveguide mode at the end of the chain does 
not show significant difference for two cases - when the 
first or the first two atoms are excited. For the rest of 
the part, we consider the case when the only first atom 
is externally excited. 

Our numerical calculation suggests that the waveguide 
mode can be achieved for a spacing a ^ 0 . 6 Ao/ 27 r. Since 
the resonance is intrinsic in the atomic system , the wave¬ 
length of resonant light Aq is species-specific. In Fig. |^we 
show the electric field magnitude, normalized by external 
electric field amplitude as a function of the position 
in the unit of Aq for different lattice spacing a. Only for 
a around 0 . 6 Ao/ 27 r is there a “waveguide” mode. For the 
other array spacing, the wave is non-propagating down 
the chain. The real part and the imaginary part of the 
field as a function of atom position are compared in the 
inset of Fig. (b). As expected, they are of compara¬ 
ble magnitude and exhibit the same spatial dependence. 
This figure illustrates the existence of a ” carrier wave” of 
high periodicity, which is usually modulated for different 
situations. 

Waveguide eondition Let be the transverse wave 
vector normal to the propagation direction. To suppress 
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FIG. 2: (color online) Waveguide mode along the linear 

chain of atoms, (a) An array of atoms is placed on a straight 
line with the separation of a. Here, the amount of position 
fluctuations is indicated by ga. (b) The normalized intensity 
of the electric field |F^| is shown as a function of the position of 
atoms for different lattice spacing a — 0.3Ao/27r (open circle), 
a = 0.6Ao/27r (solid circle), a = 0.9Ao/27r (solid triangle), and 
a = 1.2Ao/27r (solid square). Here, Eq is the the amplitude 
of the external electric held exciting the hrst few atoms. The 
inset shows the real and imaginary part of electric held for 
the propagating mode when a = 0.6Ao/27r. (c) Distribution 
of the electric held as a function of the position is shown 
for different incoming angular frequency normalized by Eq. 
(d) Fourier transform of the electric held along the progating 
direction (x) as a function of the wave vector k in units of 
27r/ (aN) for N=320 atoms on a straight line of lattice spacing 
a = 0.6Ao/27r. 


the radiation loss, k± must be imaginary. The transverse 
wave vector k± is given by k± = y/(27r/A)^ — with A 
the wavelength of the guided wave and kx the wave vector 
along propagation direction. The Fourier transform of 
the electric held in the transverse mode is shown as a 
function of the wave vector k in units of 27r/{Na) for 
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FIG. 3: (color online) (a) The electric field amplitudes are 
shown as a function of the position for N=320 atoms of dif¬ 
ferent amount of positional fluctuations: 77 = 0 (black solid 
circle)and 77 = 0.2 (red cross). In (b), the Fourier transform of 
transverse electic field is shown in units of 2'Ka/N for N=320 
atoms. 


N=320 atoms (see FigJ^ (d)). The incoming wave has 
a wave vector ki = ujNa/{2 t:c) = 35.9 in the present 
unit. This is much less than the dominant peak at /c = 
284 and therefore ky is imaginary. This large periodicity 
is mainly due to the “carrier wave”, which sustains the 
subwavelength transport. There is another peak in the 
Fourier transform at ki^ which reflects the contribution 
from atoms at the beginning of the chain. 

As kx takes its maximum at the Brillouin zone bound¬ 
ary Kb = 7 r/a for a chain of particles, the necessary 
condition to support the waveguide modes is [22j [23] 

a < A/ 2 . (3) 

Nanoparticles are required to be packed as closely as 
possible regardless of waveguides made of metallic or di¬ 
electric particles [3 Hi HU. For atoms, however, the 
resonance is intrinsic, there is much greater flexibility in 
the choice of the lattice constant. 

The above discussion shows that the one-dimensional 
lattice spacing cannot be too big. Our numerical result 
indicates that the lattice spacing cannot be arbitrarily 
small as well. Physically, as the lattice spacing becomes 
small, it becomes possible to describe the system with an 
effective dielectric constant e^f / that becomes bigger and 
bigger as the lattice constant is decreased. Since ky = 
uj‘^eeff/c^ — k‘^, it becomes more and more difficult to 
render ky negative to localize the electromagnetic wave. 
Mathematically we find that the imaginary part of the 
resonance frequency increases when the lattice spacing 
becomes very small. It becomes impossible to reduce the 
radiative loss. The details of this is described in Ref. [ 20 ] . 

Frequency and Length Dependence Now we examine 
the electric field intensity along the chain for different 
angular frequencies normalized by the resonance angular 
frequency of the atom as shown in Fig. (c). The max¬ 
imum amplitude occurs around cj/cjo = 1.175. As the 
frequency moves away, the mode becomes nonpropagat¬ 
ing, as is illustrated. We also investigated the transport 


through atoms for straight lines of two different lengths 
but the same lattice spacing. The “carrier wave” peri¬ 
odicity remains the same, but the response is bigger for 
the longer chain. Generally, the resonance energy levels 
form a discrete set for a finite chain of atoms. The spac¬ 
ing between the levels becomes smaller as the length is 
increased. It is easier to find a resonance energy level 
closer to the driving frequency. 

Position fluctuation We next analyze the possibility 
of transport through atoms with fluctuations in their po¬ 
sitions so that the deviation of the atom from its perfect 
position is equal to Here ^ is a random number be¬ 
tween -1 and 1 and thus 77 measures the amount of the 
fluctuation. The dependence of the electric field E on 
distance and its Fourier transform is shown in Fig. [^for 
different amount of randomness. The response is quite 
robust. The waveguide mode persists for all the cases we 
tested. We think this is made possible because the spac¬ 
ing of the lattice is much less than the wavelength and 
that the resonance is intrinsic and is not from a Bragg 
reflection. 

To achieve the waveguide mode along the atomic chain, 
the controlled patterning of trapped atoms is critical. It 
can be done either by locating nanoparticles lithograph¬ 
ically [25] or through the self-assembly mechanism [26] . 
In both methods, atoms can be confined within the sub¬ 
wavelength trap associated with the near-held conhne- 
ment. Another promising way is to conhne atoms in a 
hollow-core hber mm] or in a two-dimensional photonic 
crystal 13 US] producing subwavelength lattice potentials. 

In our model, we consider an array of atoms retaining a 
single two-level optical transition for supporting a waveg¬ 
uide mode. Among alkali or alkali-earth-like atoms, for 
example, ytterbium or strontium isotopes have a single 
principle optical transition due to the total electron spin 
of zero [28|. For ytterbium atoms, an ideal separation 
between atoms found above is 0.6Ay5/27r ^40nm (55nm) 
where the ^Pi Sq — ^Pi) optical transition has 

the wavelength Xyb of 399 nm (556 nm). Our analysis 
may be extended to the array of atoms with more than 
one principle resonance such as rubidium, sodium and 
potassium. 

Circular chain Finally we consider “cavities” made 
from atoms arranged on a circle. For geometric reso¬ 
nances, we have not been able to construct circular cav¬ 
ities with small radiative loss. For intrinsic resonances 
with smaller lattice spacings, this has now become possi¬ 
ble. The distribution of the electric held intensity in real 
space along the circle is shown in Fig. (c) for diherent 
frequencies with the external held along diherent direc¬ 
tions. The resonance frequency for propagation is now 
shifted from that for the straight line, the dependence 
on the resonance frequency is quite sharp. The mode is 
excited no matter what direction the external held is. 

In Fig. (b) we show the Fourier transforms of the 
electric helds in the radial and the tangential direc¬ 
tions for N=640 atoms separated by a lattice spacing 
a = 0 . 6 Ao/ 27 r on the circumferences of circle after the 
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atoms connected by a one-dimensional line waveguide [6] 
or subwavelength-diameter nanofibers m reveal inter¬ 
esting transport properties with reduced radiation. 

In conclusion, we demonstrate that the electromag¬ 
netic energy can be efficiently transported along a chain 
of atoms through an atomic resonance with a wave¬ 
length of Aq. We numerically analyze the waveguide 
mode propagating along a finite chain with appropriate 
lattice spacing a ^ 0 . 6 Ao/ 27 r. The efficient subwave¬ 
length transportation of electromagnetic energy along 
the one-dimensional waveguide would provide a main el¬ 
ement in integrated nanoscale optical circuits. We be¬ 
lieve that ultracold atoms can be localised with subwave¬ 
length lattice spacing in the plasmonic nanosphere traps 
or hollow fibers. Our study of ID array structure rep¬ 
resents a basic element of complex geometries such as 
two-dimensional patterning of atoms. For example, it 
should be interesting to study how the light propagates 
in the two-dimensional triangular array from one end to 
the other m- 

The work was supported by the Hong Kong Research 
Grants Council (Project No. ECS26300014). STC 
thanks the hospitability of the Institute for advanced 
studies and the physics department of Hong Kong Uni¬ 
versity of Science and Technology where this work was 
carried out. 


FIG. 4: (color online) Waveguide mode on the circular 

atomic chain, (a) The atom at the top is excited by an exter¬ 
nal electric field resulting propagating modes in both direc¬ 
tions. (b) Fourier transform of radial (\Er{k)\) and tangential 
{\E(f){k)\) electric fields as a function of the wave vector k in 
units of I'KajN for N=640 atoms on the circumferences of cir¬ 
cle at an angular frequency uj = 1.1. (c) Distribution of the 
radial electric field as a function of the position of the atoms 
normalized by the lattice constant a = 0.6Ao/27r for atoms 
on the circumferences of a circle with 640 atoms at different 
frequencies. 


atom at the top is excited by an external field at an an¬ 
gular frequency u; = 1.1. As can be seen , the field is 
mainly along the radial direction. We have also studied 
the system for various numbers of atoms. We show in 
Fig. El (b) the Fourier transforms of the electric fields in 
the radial direction for a smaller number of N=320 atoms 
separated by the same lattice spacing a = 0.6 on the cir¬ 
cumferences of circle after the atom at the top is excited 
by an external field. The magnitude of the field is now 
less as is expected. 

We emphasize that even though this looks similar to 
the Whispering gallery mode [29] , there is no outside wall 
to confine the wave in the present case. Indeed, there is 
a high wave-vector “carrier wave” which makes the cav¬ 
ity possible. In addition, the transport of photons with 
reduced radiative loss has been investigated also in dif¬ 
ferent systems. For example, a system of two distant 
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Supplementary Information 


2. Multiple scattering 


1. Coupling to atoms 


We have the usual Hamiltonian H = — er • E. The elec¬ 
tric field can be expanded in terms of vector spherical 
harmonics: 

= E E (4) 

Z=0 m=—l 

The vector spherical harmonics are defined as 

(5) 

Ni;ll{x,0,^) = xM^^l{x,0,cf>) (6) 

(7) 

The scalar spherical functions are given by 

0, <j>) = 4-^) (x)P4(cos (8) 


where 



J = 1 
J = 3 


(9) 


The response of the lines of atoms to an external elec¬ 
tromagnetic field is treated by multiple scattering the¬ 
ory. The total incoming electric field driving the atom 
is the sum of the external electric field E?^^ and the 
scattered fields from all the other atoms: 

E-^^ = E?^^ + This incoming field is scat- 

tered by the atom and produces a scattered field E^ ^ 
that is equal to the t matrix times the incoming field: 
E^^^^ = The field E^^^^^is equal to the photon 

propagator G^j times the scattered field from site j: 

E-^^^ = GijEj^\ To summarise, we get the set of lin¬ 
ear equations 

Ef ^ GytEj^^ = tE4* (10) 

The photon propagator [32] in terms of the vector po¬ 
tential is Because 0 = 0, we 

need not worry about /i = z/ = 0. In classical E&M, we 
worry about the Green’s function G(r) = / d^r/{k^ 

This is identical to the Fourier transform of the photon 
propagator. 


3. Small lattice spacing 


with X = kr and n = 0,1, 2,..., m = 0, ±1,..., ±n. The 
superscripts J=1 and J=3 correspond to incoming and 

outgoing waves. a;o = kro, 7mn = 

The electromagnetic field locally is a linear combination 
of M oc and oc V x These functions all 

have total (orbital plus spin) angular momentum 1. 

The Hamiltonian is = —er-E. We have r-VY/rn = 0? 
r • LYim = 0, r • rYim = r'^Yim- Thus r • Ni„ is the only 
nonzero term in the dipole approximation. M oc LYJ^ 
contains no radial component, thus this term is not cou¬ 
pled to the atom. To illustrate the physics we consider 
atoms that exhibit a p wave resonance so that the atom 
goes from a state with total angular momentum 0, |0 > 
to a state with angular momentum 1, |l,m > and then 
reemits a photon. The tangent of the corresponding scat¬ 
tering phase shift exhibit a frequency dependence given 
by tan = 1/{uj — ujq ^ iV). The expansion coefficients of 
the incoming and outgoing fields are related by the t ma¬ 
trix : = taT). Because the Hamiltonian only involves 

the electric field and not the magnetic field, there is no 
coupling involving the coefficients a^^). Since we assume 
a ground state of zero angular momentum and since the 
Hamiltonian is rotationally invariant, the t matrix is di¬ 
agonal in /,m. The t matrix is related to the tangent of 
the scattering phase shift 6 hy t = tan (tan ^ + i). For 
this resonance the coupling to the electric field is propor¬ 
tional to < 0|E-r|l,m > . Since 'Nim is of total angular 
momentum 1, by the Wigner Eckert theorem, only the 
terms with 1=1 comes in. 


In this appendix, we describe the physics when the 
lattice spacing is small. The multiple scattering equation 
is: 

e 4 = tiE ciyEf ^ + Efl (11) 

We consider a line of atoms line up along the z axis. 
Then the Green’s function is particularly simple [33]. G 
is diagonal in m and is given by 

min{n,u) 

Gy = Co E i^C;a,hMRij\) 

q=0 

for R on both the left and the right, where p = n-\-u — 2q^ 
Cp = n{n + 1) + + 1) ~ P{P + 1), a is the Wigner 3j 

symbol. For the present approximation, n = v = 1. We 
thus get {C 2 = 4 — 6 = —2, Cq = 4. 

Cy = Co[2aoh2{koRij) + iaihoikoRij)] 

, For an infinite chain, Ei ex ex.p{ikla). We need 

G{k) = Y G{R) = C'o[2ao52(fc) +4ai5o(fc)] 

R^O 

where gn{k) = '^j^^Qex.p{ikR)hn{koR) The above equa¬ 
tion becomes 


E{k) = t[G{k)E{k) -h E^^\k)] 


( 12 ) 
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The resonance condition is 

l=tG 


dx{ex.p[i{k + ko)x] + ex.p[i{—k + ko)x])/{koxa). 


We look at the small a limit. Take h = j -\- iy. 
ho{x) = (sin(x) — i cos{x))/x = —i[i sin(x) + cos{x))/x = 
—iex.p{ix)/x Define Xj = ja. For example, for n=0, 
ho{x) = —iex.p{ix)/x, 

go{k) = 2 E(A Xj) ex.p{ikxj)ho{ko\xj\)/a 

= J2{Ax j){e-xp[i{k^ko)xj]-^ex.p[i{-k^ko)xj])/{koXja) 
j>0 


Both the real and the imaginary part of g scales as 
ln{a)/a. This dependence will also dominate 1 — tG. 
Thus the imaginary part of Eg = tEij^l — tG) will not 
be small in the small a limit. 


[1] M. Bajcsy, S. Hofferberth, T. Peyronel, V. Balic, 
Q. Liang, A. S. Zibrov, V. Vuletic, and M. D. Lukin, 
Physical Review A 83, 063830 (2011). 

[2] M. Bajcsy, S. Hofferberth, V. Balic, T. Peyronel, 
M. Hafezi, A. S. Zibrov, V. Vuletic, and M. D. Lukin, 
Physical Review Letters 102 , 203902 (2009). 

[3] C. L. Hung, S. M. Meenehan, D. E. Chang, O. Painter, 
and H. J. Kimble, New Journal of Physics 15, 083026 
(2013). 

[4] A. Goban, C. L. Hung, S. P. Yu, J. D. Hood, J. A. Muniz, 
J. H. Lee, M. J. Martin, A. C. McClung, K. S. Choi, D. E. 
Chang, et al.. Nature Communications 5 (2014). 

[5] D. Dzsotjan, A. S. S 0 rensen, and M. Eleischhauer, Phys¬ 
ical Review B 82, 075427 (2010). 

[ 6 ] J. T. S. S. Pan, Optics Letters 30, 2001 (2005). 

[7] D. E. Chang, L. Jiang, A. V. Gorshkov, and H. J. Kimble, 
New Journal of Physics 14, 063003 (2012). 

[ 8 ] E. Ozbay, Science 311, 189 (2006). 

[9] W. L. Barnes, A. Dereux, and T. W. Ebbesen, Nature 
424, 824 (2003). 

[10] R. Zia, J. A. Schuller, A. Chandran, and M. L. 
Brongersma, Materials today 9, 20 (2006). 

[11] S. A. Maier, P. G. Kik, H. A. Atwater, S. Meltzer, 
E. Harel, B. E. Koel, and A. A. G. Requicha, Nature 
Materials 2, 229 (2003). 

[12] S. 1. Bozhevolnyi, J. Erland, K. Leosson, P. M. W. Skov- 
gaard, and J. M. Hvam, Physical Review Letters 86 , 3008 
( 2001 ). 

[13] P. Berini, Physical Review B 61, 10484 (2000). 

[14] J. R. Krenn and J. C. Weeber, Philos.Trans.R. 
Soc.London, Ser.A 362, 739 (2004). 

[15] K. Tanaka and M. Tanaka, Applied Physics Letters 82, 
1158 (2003). 

[16] S. A. Maier, Plasmonics: fundamentals and applications: 
fundamentals and applications (2007). 

[17] J. Du, S. Liu, Z. Lin, J. Zi, and S. T. Chui, Physical 
Review A 79, 051801 (2009). 


[18] M. Gullans, T. G. Tiecke, D. E. Chang, and J. Eeist, 
Physical Review (2012). 

[19] A. Gonzalez-Tudela, C. L. Hung, D. E. Chang, J. Cirac, 
and H. J. Kimble, Nature Photonics 9, 320 (2015). 

[20] See supplemental material (????). 

[21] P. C. Waterman, Proceedings of the IEEE 53, 805 (1965). 

[ 22 ] A. L. Burin, H. Gao, G. C. Schatz, and M. A. Ratner, 
JOSA B 21 , 121 (2004). 

[23] G. S. Blaustein, M. 1. Gozman, O. Samoylova, 1. Y. 
Polishchuk, and A. L. Burin, Optics Express 15, 17380 
(2007). 

[24] M. Quinten, A. Leitner, J. R. Krenn, and E. R. 
Aussenegg, Optics Letters 23, 1331 (1998). 

[25] P. Nagpal, N. C. Lindquist, S.-H. Oh, and D. J. Norris, 
Science 325, 594 (2009). 

[26] J. A. Pan, C. Wu, K. Bao, J. Bao, R. Bardhan, N. J. 
Halas, V. N. Manoharan, P. Nordlander, G. Shvets, and 
E. Capasso, Science 328, 1135 (2010). 

[27] C. A. Christensen, S. Will, M. Saba, G.-B. Jo, Y.-i. Shin, 
W. Ketterle, and D. Pritchard, Physical Review A 78, 
033429 (2008). 

[28] J. E. Sansonetti and W. C. Martin, Journal of Physical 
and Chemical Reference Data 34, 1559 (2005). 

[29] L. Rayleigh, Rayleigh: The Theory of Sound, 1894 
- Google Scholar (Republished by Dover Publications, 
1945). 

[30] L. Tong, R. R. Gattass, J. B. Ashcom, S. He, J. Lou, 
M. Shen, 1. Maxwell, and E. Mazur, Nature 426, 816 
(2003). 

[31] C. Becker, P. Soltan-Panahi, J. Kronjager, S. Dorscher, 
K. Bongs, and K. Sengstock, New Journal of Physics 12 , 
065025 (2010). 

[32] E. Mandl and G. P. Shaw, Quantum Field Theory (John 
Wiley V Sons, 2013). 

[33] Y.-l. Xu, Applied Optics 34, 4573 (1995). 



